Two parameters scaling approach to Anderson localization of weekly interacting BEC 
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We numerically study the Anderson localization of weekly interacting Bose-Einstein condensate in a one- 
dimensional disordered potential. We show that two parameters are needed to completely describe such system, 
and the density profile of which can be described with the sum of two exponential functions. This is a new 
attempt for precise description of systems with interplay of disorder and interaction. 
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I. INTRODUCTION 

Disorder is ubiquitous in nature that strongly affects the 
properties of many physical systems even it is only a weak 
perturbation. Fifty years ago, the localization of individual 
particles or waves in a disordered crystal was predicted by An- 
derson ill], and thus it is called Anderson localization (AL). 
It can be understood as an effect of multiple refection of a 
plane wave subject to random scattering or random potential 
barriers 101 ■ Later, AL has been unambiguously observed in 
many systems under the approximation of a single particle in a 
stationary disordered potential. For exarnple, electromagnetic 
waves in photonic lattices with disorder olQl and momentum 
distribution of quantum kicked rotor 15|-Ial- In single-particle 
approximation, the interaction among individual particles is 
not taken into consideration. However, real materials go far- 
beyond this approximation ifiol [Till , and thus observing AL 
is difficult. This can be understand as follows. Go beyond 
single-particle approximation, one needs to consider the in- 
teraction among particles. When this interaction induced non- 
linear effect presences, the reflected waves will interfere with 
each other. Therefore, fully understanding the interplay of 
disorder and interaction is an extremely difficult task both ex- 
perimentally and theoretically 1 1 ij. 

Ultracold quantum gas possesses unprecedented possibility 
of controlling almost all relevant physical parameters 1I12IJ17I1 . 
and thus recognized as an ideal system for quantum simula- 
tion. A particularly interesting aspect of the system is im- 
plementing random speckle potential using laser beams lisl - 
I22II . which makes this system suitable for observing AL. Re- 
cently, take advantage of highly controllability of this sys- 
tem, two experimental groups have reported the observation 
of localization of a noninteracting as well as weak interact- 
ing Bose-Einstein condensate (BEC) in two different kinds of 
disordered potentials [23, 241. The final state profile is the- 
oretically described by a single parameter called localization 
length (LL) ll26lj30ll . However, when a BEC stop to expand 
in a disorder potential, its interaction energy is not completely 
converted into kinetic energy. Therefore, one can successfully 
describe the wing of the localized profile in terms of LL under 
the single-particle approximation as there only has negligi- 



ble particle interaction due to long time free expansion. But 
this does not applicable for the center of the localized profile 
where interaction between particles can not be ignored. This 
fact naturally leads us to consider the possibility of consider- 
ing the problem of describe AL in interacting system with two 
parameters, one for the wing and the other for the center. 

In this paper, we consider a concrete example of one- 
dimensional (ID) BEC with repulsive interaction in a random 
potential. We prove that two parameters description of AL in 
such system is more reasonable than that of single LL param- 
eter. This is a new method to describe the localized profile, 
which provides a new method for studying AL with the in- 
terplay of disorder and interaction. The rest of this paper is 
arranged as following. In section II, we describes ID BEC 
system with repulsive interaction, and shows that there is two 
different LL for the wing and center parts. In section III, as 
the wing LL is well known, we present detail study of the cen- 
ter LL focusing on its scaling law. In section IV, we give an 
approximate analytic expression linking the density profile of 
atom to the two LL and discuss the deviation of the expres- 
sion. Finally, a brief summary is given in section V. 



II. TWO PARAMETERS DESCRIPTION OF AL 

Let us consider ID BEC with repulsive interaction initially 
loaded in a harmonic trapping potential Vhoiz) - jMoip-z', 
where M is the atomic mass and o) is the trapping frequency. 
The effective ID structure can be achieved by applying an 
extremely tight harmonic vertical confinement to froze the 
atomic motion in the other two dimensions. The atomic in- 
teraction is effectively characterized by the i-wave scattering 
with effective ID interacting strength labeled as g, which is 
experimentally tunable using the Feshbach resonance tech- 
nique Isill . Here, we consider BEC in weakly interacting 
regime, i.e. n » Mg/h^, where n is the average atomic den- 
sity. Under the mean-field approximation, the dynamics of 
the considered system is governed by the following Gross- 
Pitaevskii (GP) equation 
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where pz is the momentum operator, ju is the chemical poten- 
tial, Vext(z) is the whole external potential, and the normalized 
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wave function \jj{z, t) is corresponding to a constant total num- 
ber of atoms N = J \i//iz, t)\~dz. 

Now we look into the time-evolution of the BEC in a disor- 
dered external potential governed by Eq. ([TJ. To this end, we 
follow the procedures in the experiment 1,24,] : Firstly, prepare 
the ID BEC at equilibrium of the harmonic potential with- 
out disorder, i.e., Vg^,(z) = V/uAz)- In the Thomas-Fermi (TF) 
regime (pi » hai), the initial wave function takes the form of 
an inverted parabola | fl3l l 
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(2) 



where Ljf - yj2ijj(ji? is the TF half-length, and denotes 
the Heaviside step function. Secondly, at time f = 0, one 
switches off the harmonic potential and applies a disorder po- 
tential along the expanding axis (i.e. z axis) of the BEC. We 
also assume that the disorder potential Vdiz) is generated by 
laser speckle method as in the experiment p^fl. It is a ran- 
dom potential with a truncated negative exponential single- 
point distribution ll25ll : 



Vr \ Vr 



(3) 



The average of disorder potential is set to be {Vd{z)) — 0, and 
its con-elation function C{z) - {Vd{z')Vd{z' + z)) - V^c(z/cr^) 

with c(m) - sir?{u)/u^ and Vr - -^JiV^} being the standard 

deviation. Here C{z) characterizes the amplitude of the disor- 
der fluctuation and the correlation length of the disorder Thus 
the external potential after f = is given by Ve^,{z) - Vdiz). 
In the case of o-r/^i,,, < 1 with ^,„, - -y/4Myu being the initial 
healing length of BECs, the density profile of the BEC will 
take a form of an exponential-decay function when evolving 
enough time ll24ll . 

To investigate the AL of the system, we work out the time 
evolution of the atomic density profile n{z, f) - \^{z, f)l^ with 



ipiz, f) -T exp 



'GpdtWz,Q), 



(4) 



where the GP Hamiltonian is Hcp - ji^ + Vd{z)+g\ilr{z, t)\~ and 
T being the time ordering operator We numerically calculate 
f) b y us ing the standard operator- split method. According 
to Ref. 1I32I1 . Eq. (|4]i can be rewritten as 



X exp J- J + 0(6P)^ i/fiz, t). 



(5) 



where the high-order term 0(3P) comes from the non- 
commuting relation of the terms in Hcp. In the sufficiently 
short time step 6t, this term can be safely neglected. Combin- 
ing with the Fourier transform between the position and mo- 
mentum spaces, we can finally get the numerical solution of 
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FIG. 1: (Color online) (a) The stationarity of the localized pro- 
file with disorder potential strength Vr = 0.15//, nonlinear inten- 
sity Ng/h = 0.035 Hz, total number of atoms A' = 1.7 x 10'' and 
CrI'^m ~ 0.65. There are clearly two different LLs. (b) Central LL 
versus amplitude of the disordered potential in different nonlinear in- 
tensity: A'^//!=0.035Hz (blue x), 0.104Hz (red V), 0.204Hz (green 
o). Insert: Wing LL versus amplitude of the disordered potential in 
different nonlinear intensity: A'g//?=0.035Hz(blue o), 0.104Hz (red 
□), 0.204Hz (green «). Both LLs are obtained in an evolution time 
t = 2s. 



^{z, t) following the computation procedure step by step with 
time step St. 

As expected, as shown in Fig. 1(a), the final density profile 
of BEC takes the form of similar exponential function when 
evolving in disorder potential. We find that the wing of density 
profile can be exponentially fit by a wing LL denoted as Lw in 
Fig. 1(a), which has been considered as the single parameter 
to characterize LL of the system in experiment ll24ll and takes 
the analytic form of |[l2ll 



7:M^VlcTR^l,(l-(TR/^i„,y 



(6) 



However, both in our simulations and in experiment ll24ll . the 
central part of the density profile can not be well-characterized 
by Lw though it is also localized. To describe the AL of the 
central part of the BEC, we introduce a center LL denoted as 
Lc as shown in Fig. 1(a). Here Lc is defined by the length 
in z axis which reduces the maximum value of the finial wave 
function by a factor of 1/e. Furthermore, by varying system 
parameters, such as disorder strength, interacting strength, 
etc., Lw + Lc is always hold in our numerical simulation, 
as shown in Fig. 1(b). In particular, we study the effects of 
the nonlinear atomic interaction on the AL of the two parts, 
i.e. the Lw and the Lc, in varying disorder potentials. To 
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this end, we look into the trend of two LLs by changing the 
ampUtude of the disordered potentials in different nonlinear 
interaction intensity Ng and by keeping other parameters un- 
changed. The results, as shown in Fig. 1(b), shows that the two 
LLs both nearly exponentially decrease with the increasing of 
the amplitude of the disordered potentials, but interestingly, 
they exhibit different behavior: Lw is insensitive to the non- 
linear intensity as expected from Eq. ^ while Lc is strongly 
affected by it. This can be understand by the facts: i) for larger 
amplitude of the disordered potential AL is more significant, 
and thus LLs are both smaller; and ii) the impact of g\i//(z)\~ 
on the wing part Lw is slight, but they can be significant for 
the center part Lc since the density of this part is much larger 
than that of the wing part (cf. Fig. 1(a)). 

Concluding this section, we have shown that it is better to 
use double LLs, the traditional wing one Lw and a newly de- 
fined center one Lc, to completely characterize the AL proper- 
ties of the ID weakly interacting BECs in a disorder potential. 
This double LLs description is quite necessary if we consider 
the effects of atomic interaction on the AL of the system since 
the two LLs have different responses on the interaction. 
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ni. A SCALING LAW OF THE CENTER LOCALIZATION 
LENGTH 

In contrast to the single-parameter scaling theory for the 
disorder system in finer regime, fiere we need two LLs as 
localization parameters to completely describe this nonlinear 
disorder system. The key in this nonlinear system is that the 
BEC finally reach nearly equilibrium in the disorder potential, 
and at that time the wing part will fully expand and then it can 
be well described with the LL for the noninteracting cases; 
however, the central part still contains certain residual inter- 
action energy without fully expanding, and thus the LL of this 
part could be different. 

To further investigate localization properties of the center 
part of the BEC, we numerically calculate Lc as a function of 
the nonlinear interaction intensity Ng with disorder potentials 
fixed. As shown in Fig. 2(a), we find that Lc approximately 
linearly depends on '\JNg, i.e. 



(7) 



Lc + constant. 



for a disorder potential with fixed Vr. For small Vr cases, 
there are some variances in the fitting (cf. Fig. 2(a)). Consid- 
ering the particle number with normalization equation 



\ifr{z,Ofdz = 

OO J-Ljf 8 



1 - 



[2 



3^ 



(8) 



we can obtain the expression of the TF half-length as Ltf - 
3Ng/4-fi, which implies that one has the relation Lc oz yfLjF. 
Combining with the dimensional analysis, it is natural to guess 
that there may be a relationship between Lc and ^|LTF x ^iLw- 
Through a simple fitting, we find that Lc has an interesting 
scaling law of 



FIG. 2: (Color online) (a) The Lc versus in different V^: Vi;=0. lu 
(blue *), 0.2u(red +), 0.3u(green x) with black lines are theoretical 
predicted results by Eq. (6). (b) The Lc versus the \ILv/Ltf- They 
are obtained different nonlinear intensity: Ng/h=0.035Hz (blue x), 
0. 104Hz (red *), 0.204Hz (green +), and different evolution time: 
2.1s (gray o), 4.5s (purple □), 7.0s (dark green V). The solid line is 
the function v = x. 



as shown in Fig. 2(b). 

Here we do some analysis about the physics picture of the 
scaling law (|9]l. Lw can be taken as a main measurement of 
the strength of localization of the system, thus Lc as an addi- 
tional measurement should be positively related to Lw- Since 
the central part of the BEC can not fully expand and has less 
kinetic energy, it is less sensitive to the strength of localiza- 
tion than the wing part does. Therefore, its power index is 
less than 1 and exhibits Lc ^ yfLw- Additionally, Ltf is a 
measurement of the size and the profile of the initial state, 
and then Lc could also be positively related to Ljf- Because 
of the disorder potential and the smaller kinetic energy, the 
power index is also less than 1 and exhibits Lc <^ y/LjF- 



IV. AN APPROXIMATE ANALYTIC EXPRESSION OF 
DENSITY PROFILE 

Now we turn to look for an analytic expression of density 
profile of the localized BEC using the two LLs. The previous 
analysis shows that its finial density profile takes the form of 



n(z, t) ' 



exp(-2|z| /Lc) kl < ZQ 
exp(-2|z|/Liy) |z| > zq 



(10) 



Lc - yjLwLj 



(9) 



where t is the typical time scale when the BEC is nearly sta- 
ble, and zo is the theoretical cross-point. Here zo is diflicult to 
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FIG. 3: (Color online) The scale parameter L^/Lrf versus the resid- 
ual interaction energy AE with iTr/§;„, = 0.65. Lw/Ltf is obtained 
from different nonlinear intensity: Ng/h=0.035Hz (blue ♦), 0.104Hz 
(red *), 0.204Hz (green +), and different evolution time: 2.1s (gray 
o), 4.5s (purple □), 7.0s (dark green V). The solid line is the best fit 
with the function in the Eq.(lO). 




FIG. 4: (Color online) The parameter 6z versus the residual interac- 
tion energy AE and Zi/Lc- The former is the solution of Eq. l ll5t . 



Substituting Eq. ( fT3] l into Eq. (fT2l i. one can work out A, 
which is approximately given by 



determine, so we assume a simple form 

naiz, t) = exp(-2\z\/Lw) + exp(-2k|/Lc) (1 1) 

to describe the full density profile, where A is an undeter- 
mined coefficient. In the following we will see that A can 
be approximately determined. 

Let us consider the normalized residual interaction energy 
AE, which is defined as 



A£ 



g j mz,OWz 

5 
6A^ 



, , , 4(A^-A) 
+ (N- A) Y + 



(12) 



A ^ 



0.75 1 ^\ + 0.04 



N 



(14) 



Up to this, we have shown that the full density profile of the 
localized BEC can be characterized by an approximate ana- 
lytic function Uaiz, t) as shown in Eq. (fTTl i. where the two 
LLs are used and the coefficient A can be determined by the 
equation ( fT4] i. 

Now we investigate the deviation in the determining Lc by 
using Eq. ( fTOb and Eq. ( fTTI ). To this end, we define z\ and zi 
{z\, Z2< Zq) as the solutions of n(zi, t) = riaizi, t) in the central 
part (i.e. z\,2 < Lc), such that we have 



N-A A 

H 

Lw Lc 



-2-p- 



N-A _2ii A _2iJ- 

e '■»' H e '-<: . 

Lw Lc 



(15) 



with p - ^ILtfILw- Note that we have adopted the scaling 
law (|9]) in the derivation of Eq. (fT2l i. For the single LL cases, 
i.e. A = OinEq. (fT2li. we can find AE - 5Ljf/6Lw- Based on 
the numerical simulation, we also find that there is a well fit- 
ting formula to characterize the relationship between AE and 
L^./Ltf for our two LLs description, as shown in Fig. 3. The 
fitting formula is given by 



Ltf 

This relation can be understood by the fact that Lw and Lj-f 
are the characteristic lengths of finial and initial states, and 
thus their ratio may has certain connection with the residual 
interaction energy which is the ratio of the interaction energy 
of final and initial states. If the BEC expands fully, the resid- 
ual interaction energy tends to zero. Then, Lw will be much 
larger than Ltf, corresponding to AE for Lw/Ltf °° 
in Eq. (fTSl l (cf. Fig. 3). On the other hand, if the expan- 
sion is relatively very small, comparable to the initial length 
scale of the BEC, then Lw will be nearly equivalent to Ltf, 
corresponding to the residual interaction energy tends to unit. 



We then define the deviation 6z as 

(16) 

Zl 

& > indicate the distribution of our approximation is wider 
than the precise distribution in the z-direction , and vice versa. 
According to Eqs. ( fT3l fT4l [T5] ). we numerically compute the 
deviation 5z as a function of AE and z\ ■ The result in Fig. 4 
shows that the deviation of z-direction with the same n{z, t) 
and parameters. And it can find that when AE is in the small 
and large sides the deviation is smaller than that in the inter- 
mediate regime. However, the deviation using our approxima- 
tion is always less than 18% in the whole regime. Relatively 
speaking, for non-interacting system, both the deviation of the 
experimental and theoretical values are larger than 50% ll24ll . 
While for the interacting system, the deviation predicted in 
Ref. lfl2ll is larger than 60% with the same parameters of Fig. 
(1) when Vr is in the small and large sides. Therefore, we 
can conclude that Eq. ( fTTT i of is a better approximation to 
describe the density profile of the localized BEC in a wide 
parameter range. 
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V. CONCLUSION 
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In summary, we have demonstrated that it is better and more 
completely to use double LLs to describe the AL of ID weekly 
interacting BECs in a disordered potential. We furthermore 
find a scaling law related to the relationship between the newly 
defined LL and the nonlinear atomic interactions. An approx- 
imate analytic form of the full density profile of the localized 
BEC is also proposed by using the two LLs. 
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